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We present a new astro-physical site of the big bang nucleosynthesis (BBN) that are very peculiar 
compared with the standard BBN. Some models of the baryogenesis suggest that very high baryon 
density regions were formed in the early universe. On the other hand, recent observations suggest 
that heavy elements already exist in high red-shifts and the origin of these elements become a big 
puzzle. Motivated by these, we investigate BBN in very high baryon density regions. BBN proceeds 
in proton rich environment , which is known to be the p-process like. However, by taking very 
heavy nuclei into account, we find that BBN proceeds through both the p-process and the r-process 
simultaneously. P-nuclei such as 92 Mo, 94 Mo, 96 Ru, 98 Ru whose origin is not known well are also 
synthesized. 

PACS numbers: 26.35.+C, 98.80.Ft, 13.60.Rj 



I. INTRODUCTION 

What happened in the early universe has a great influ- 
ence on the history of the universe because they deter- 
mined the initial conditions. It is very important to check 
whether our standard model of cosmology is correct or 
not as theories and observations develop. Baryogenesis 
and BBN should be checked because they determine the 
history of the chemical evolution. 

In the standard model of elementary particle physics, 
the baryogenesis is possible only through the electro- 
weak spharelon process. In the supersymmctric standard 
model, it is much easier to explain the baryon number 
asymmetry because there are many scalar fields which 
have baryon number. One of the most striking property 
of supersymmetric theories is that they have flat direc- 
tions in potentials. Some of them have baryon number 
and if fields condensate in these directions, it is possi- 
ble to produce large baryon number. This is the basic 
idea of the Affleck-Dine baryogenesis [lj . Usually baryon 
number production is assumed to be taken place homo- 
geneously all over the space. This is natural because we 
know that the universe is homogeneous and if baryogen- 
esis is inhomogeneous in large scale it contradicts obser- 
vations 2] . Of course resolution ability of observations is 
limited and small scale inhomogencity is not excluded by 
observations. Though it seems to be unnatural to con- 
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sider such small scale inhomogeneity, recent observations 
force us to reconsider the possibility of inhomogeneous 
baryogenesis. 

It has become clear that the evolution of matter started 
earlier than we have known before. For example, Wilkin- 
son Microwave Anisotropy Probe (WMAP) data suggests 
that reionization began when z ~20 3]. According to 
Refs. Hill, 

star formation activity started when z > 10. 
In addition, it is known that the quasar metallicity did 
not significantly change from the time of high rcdshift 
to the present time |g. Recently a galaxy at z=10.0 was 
observed ■ Other evidences of heavy elements from the 
high redshifts are given in H H ES 111 E El 

Motivated by these observational evidences, we inves- 
tigate the possibility that inhomogeneous baryogenesis 
produced very high baryon density in small fraction of 
the universe and in these regions some fraction of heavy 
elements were already synthesized during BBN. 

Heavy elements production during BBN itself is not 
a new idea. Previous researches on the inhomogeneous 
big bang nucleosynthesis are given in fl4L[l5l . Heavy ele- 
ments production is also mentioned in |l6j |. These works, 
however, do not include very heavy elements |14| or they 
create neutron rich regions and calculate the nucleosyn- 
thesis in those regions EH- These are not suitable for 
our present purpose. 

For the production of heavy elements during BBN, 
high baryon density is necessary. However if we simply 
increase the baryon-photon ratio all over the universe ho- 
mogeneously, it would apparently contradict the observed 
light element abundances E3 and CMBR 2]. Instead, 
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we assume that the baryon density of the universe is in- 
homogeneous before and during BBN. In most part of the 
universe 77 is small (77 ~ 6x 10~ 10 ) as observed while small 
fraction of the universe is occupied with very high baryon 
density, r\ ~ 0(1). Because our aim is to see how BBN 
goes in the high baryon density regions and not make the 
precise adjustment between BBN and CMBR, we neglect 
the baryon diffusion. In this case, the baryon density in 
high density regions can be treated almost free parameter 
without contradicting observations. (It is a complicated 
problem whether we can treat 77 as a free parameter in 
realistic models. See, for example, [T^.l 

In section [H] we explain the theoretical aspects of our 
model |Tfl| . In section ITTT1 we explain our network and 
what kind of effects we take into account. Section IIVI 
is the main results of our numerical study. In this sec- 
tion we explain BBN is the p-process like and simulta- 
neously the r-process like. And also BBN can produce 
very heavy elements including proton rich nuclei such as 
92 Mo, 94 Mo, 96 Ru, and 98 Ru. 



II. THEORETICAL BACKGROUND 

Theoretical back grou nd of this model is inhomoge- 
neous baryogenesis [Tj|. We are going to explain basic 
aspects of this model. For more detail, see [Ta.l20l| . 

The basic idea of the model ^3 is a modified version of 
the Affleck-Dine baryogenesis [j| . Assume that the inter- 
action Lagrangian has the general renormalizable form 

£ mt = AH 2 $ 2 +.g|0| 2 $ 

= A($-a> 1 ) 2 |0| 2 -A$ 2 |0| 2 , 

where <j> is Affleck-Dine (AD) field, $ is the inflaton field, 
g and A are the coupling constants and <!>i = —g/2X. 

In a simplest case, the effective mass of the AD field <f> 
can be written by 

(mf // ) 2 = m 2 + A(<I>-$ 1 ) 2 , (2) 

where m 2 , is the vacuum mass of <f>. 

The vacuum expectation value of $ is assumed to 
evolve from very large value, i.e. , $ > $1, decreases 
to zero. As $ goes down to ~ $1, the effective mass 
square of cj> becomes negative and the phase transition 
takes place. When $ is far from $1 the mass square is 
positive. If the duration of $ ~ $1 is short, the transition 
would take place only in a small fraction of space. Con- 
sequently in the dominant part of the universe baryon 
asymmetry is small as observed 77 — 0(1O~ 9 ), while in a 
small part of the universe the baryon asymmetry can be 
very large, even close to unity. In this simple model, the 
signature of barionic charge is not fixed. Baryonic chaege 
can become both positive and negative [21j . However, 
the high density regions are very small compared to cos- 
mological scale, high density anti-matter regions would 



disappear by pair-annihilation while late time inflation 
can prevent the annihilation 22]. 

Because we are not very interested in the detail of the 
shape of the bubbles and the effect of diffusion in this 
paper, we assume that the bubble sizes are large enough 
to neglect the diffusion effects. Also, bubbles are not 
large so as to contradict the observations Q. 

In this case, the BBN calculation can be treated as that 
of homogeneous big bang nucleosynthesis. In the follow- 
ing section, we present the results of the calculations and 
their physical interpretations. 

At first sight, in the high baryon density regions the 
reaction seems to proceed along the proton rich side be- 
cause BBN occurs in proton rich environment |2£j • How- 
ever, surprisingly it is not correct. BBN proceeds along 
the proton rich and the neutron rich side. 

III. NUMERICAL CALCULATIONS 

The basic method of our calculation is the same as that 
of the homogeneous big bang nucleosynthesis. 
We solve the Friedmann equation 

(a\ 2 _ 8irGp 

[a) - 3 (3) 

where p = p-y + (p e - + p e + ) + p v + pb, and a is the scale 
factor. The energy conservation law is 

for the time evolution of the temperature and the baryon 
density. 

Abundance change in the region is evaluated with a nu- 
clear reaction network, which includes 4463 nuclei from 
neutron, proton to Americium (Z = 95, A — 292). Nu- 
clear data, such as reaction rates, nuclear masses, and 
partition functions, are same as in |23|. It should be em- 
phasized that both proton-rich and neutron-rich nuclei 
are produced in a high density region (Fig.^J. Therefore 
it is required using a large network to calculate abun- 
dances in the high 77 region. 

IV. RESULTS 

We have calculated BBN for various values of 77, from 
10~ 10 to 10~ 2 . It is known that in the standard, low 
baryon density BBN, nuclei heavier than Boron are 
hardly synthesized. Fig. ^ represents the synthesized nu- 
clei for 77 = 1 x KT 6 at the epoch T = 1 x 10 7 K. We 
can see that heavier nuclei such as Ca (10~ 14 in mass 
fraction) are synthesized. 

Naturally as 77 becomes large, heavier nuclei are syn- 
thesized. However we find there is drastic change in the 
nucleosynthesis around »)~3x 10~ 4 . To see this transi- 
tion, we pick up two values of 77 = 10~ 4 and 1CP 3 , and 
investigate what is going on during the nucleosynthesis. 
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(N=20) Neutron 

Fig. 1: Abundance distribution at T = 1 x 10 7 K (r) = 10" 6 ) is 
plotted on the nuclear chart. The gray color regions represent 
stable nuclei. The red color represents more synthesized nuclei 
and blue color represents less synthesized nuclei. Those nuclei 
that are not synthesized in the standard BBN scenario, such 
as Ca are synthesized. 



Fig. [2] represents how many nuclei are synthesized at 
the temperature T = 3x 10 9 K, and r\ = 1CT 4 . 

Red (blue) color represents more (less) synthesized nu- 
clei. (Stable nuclei also plotted in Fig- El 01 and [5] with 
gray color.) We can see that the reaction goes along sta- 
ble line. It is well known that nuclei whose neutron and 
proton numbers are special values (for example 20, 28, 
50, 82 etc.), the magic numbers, are especially stable. At 
these points, the reactions are stagnated and the reac- 
tion paths are bent. Especially, the stagnation at N = 
82 is one of the biggest factors that prevent the reaction 
to proceed further beyond the mass number 190. As the 
temperature goes down, the locus of the reaction begins 
to bent to different directions (Fig. 03). 

For lighter nuclei (mass number A < 100), proton cap- 
tures are very active and the locus moves to proton rich 
direction. For nuclei whose mass numbers are between 
100 and 120, the locus is across the stable nuclei from 
proton rich side to neutron rich side. For heavier nu- 
clei for A > 120, neutron capture is more efficient. This 
suggests that both the r-process and p-process occur si- 
multaneously in BBN. 

Physical interpretation of this situation is as follows. 
The environment in BBN is proton rich, i.e., the electron 
fraction Yc ranges from 0.8 to 0.9 |2jj. Naive expecta- 
tion of BBN is the p-process. For relatively light heavy 
nuclei, proton capture is active. However proton cap- 
ture processes become exponentially difficult as the pro- 
ton number increases because of their coulomb barrier. 




(Z=50) 



Fig. 2: Abundance distribution at T = 3 x 10 9 K (r) = 10" 4 ) 
are plotted on the nuclear chart. The red color represents 
more synthesized nuclei. Produced nuclei distribute along 
the stable nuclei. 



On the other hand, neutrons are still not consumed out 
during heavy nuclei are synthesized as shown in Fig. ^ 
Very heavy nuclei captures neutrons and the locus of the 
reaction changes toward the neutron rich side. Transi- 
tion point from proton rich side to neutron side depends 
on the baryon-photon ratio r\. The transition occurs at 
larger mass number for larger r\. The reactions depend 
on the abundances of the seed nuclei. The higher baryon 
density follows many seeds which lead to proton captures 
on heavier nuclei. Fig.0]shows the time evolution of mass 
fraction for nuclei whose mass number is 90 and 158. Let 
us see the time evolution of the mass number 90 . 90 Zr is a 
stable nucleus and 90 Mo is a proton rich one. First, 90 Zr 
is synthesized and later 90 Mo is synthesized while the 
amount of 90 Zr decreases. This represents that the reac- 
tion first proceeds along the stable line and later move 
to the proton rich side. At the late stage, 90 Zr increases 
while 90 Mo decreases because unstable proton rich nuclei 
decay to stable nuclei such as 90 Mo — ► 90 Nb — > 90 Zr. The 
lifetimes of 90 Nb and 90 Mo are 14.6 h = 5.33 x 10 4 sec 
and 5.67 h = 2.04 x 10 4 sec, respectively. 

For heavier nuclei of the mass number 158, the situa- 
tion is different. 158 Gd is a stable nucleus and 158 Eu is 
a neutron rich nucleus, instead of proton rich one. First 
the stable nucleus 158 Gd is synthesized and later neutron 
rich 158 Eu is synthesized. This shows that in heavier nu- 
clei region, the stable nuclei are produced first as 90 Zr, 
but later the neutron rich nuclei are produced instead 
of proton rich nuclei. The decrease in 158 Eu in the late 
stage is the same as 90 Mo, (3 decay to stable nuclei. The 
abundances of neutron rich nuclei of the mass number 
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Fig. 3: Abundance distribution at T = 1 X 10 9 K (?) = 10" 4 ). 
For nuclei with A < 100 the synthesis is the p-process while 
for heavier nuclei the synthesis is the r-process. 
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Fig. 4: Time evolution of mass fraction (fj=lx 10~ 4 ). Neu- 
trons are still left when heavy elements are synthesized. For 
nuclei with mass number 90, stable ones are synthesized first 
followed by proton rich nuclei. On the other hand, for nu- 
clei with mass number 158 stable nuclei are synthesized first 
followed by r-rich nuclei. 



90 and those of proton rich nuclei of the mass number 
158 are very small and not drawn in this figure. We can 
also see that neutrons are still left when heavy nuclei are 
synthesized. 

Now let us see the case r\ =10~ 3 . Fig. El shows the 
locus of the reaction at the temperature T = 1.8 X 10 9 
K. It is apparently different from the results of rj= 10~ 4 . 
For, in this case, the reactions first proceeds along the 
stable line. However, the reactions directly proceeds to 
the proton rich region. Another important difference is 
that very heavy nuclei of A > 80 are not synthesized. 
The physical interpretation is as follows. In a high baryon 
density region, the seeds for the reactions to proceed are 
abundant. The nuclear reaction proceeds promptly and 
all neutrons are consumed by light nuclei as shown in 
Fig- El This prevents the nucleosynthesis from proceeding 
to the large mass number region. In Fig. EJwe only draw 
the abundance having A = 90. Heavier nuclei are not 
synthesized enough. When heavy nuclei are synthesized, 
neutrons are almost consumed out. 




Fig. 5: Abundance distribution at T = 1.8 x 10 9 K (77 = 10" 3 ). 
Nucleosynthesis occurs through the p-process and very heavy 
nuclei are not synthesized. 



In Fig. [7J we show the relation between the mass num- 
ber and the number fraction relative to the solar abun- 
dances. As the baryon density becomes larger, the heav- 
ier nuclei are synthesized for r\ less than 1 x 10~ 4 . How- 
ever, when 77 > 1 x 10 -4 , the maximum mass number 
decreases as 77 becomes larger. 

The number fraction ratios of p-nuclei to the solar 
abundances is listed in Table. [I] 

For 77 = 10" 3 , 92 Mo, 94 Mo, 96 Ru and 98 Ru drastically 
increase, due to the change of the loci of the reaction 
flows. This suggests that highly inhomogeneous BBN 
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Fig. 6: Time evolution of mass fraction (77 = 1 x 10 -3 ). Neu- 
trons have already been consumed when the abundances of 
heavy elements are increasing. 
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Table. I: The number fraction of p-nuclei relative to the solar 
abundances. 

would have a large influence on the abundances of solar 
p- nuclei. 

The observed abundances seem not to be explained by 
BBN with only a single 77. However this does not ex- 
clude the possibility of our assumption. It is unnecessary 
for BBN abundances to match exactly to the solar abun- 
dances because produced nuclei in high r\ regions would 
have mixed with nuclei synthesized in low density regions 
and also there should be nuclei synthesized in star activ- 
ities. 

Basic feature of BBN at each 77 are classified as follows. 
For 77 < 1CP 6 , synthesized nuclei are limitted to A < 40. 
Nuclei whose A are around 20 are less synthesized even at 
large value of 77. For rj from 10 -5 to 10 -4 , the abundances 
of nuclei whose mass number of 30 < A < 56 grow rapidly 
with A. Abundances of nuclei A > 56 suddenly decrease 
but again slowly increase. At around A = 140, they turn 
to decrease. 

After the rapid decrease in the abundances of A > 56 
for -q = 10 -3 and 10~ 2 , the abundance profiles are rather 
different. 

For 77 = 10~ 3 , the abundances do not drastically 
change from A = 64 to around 86. They rapidly de- 
crease for A above 100 and maximum A synthesized is 
114. For -q = 10~ 2 , right side wing of Fe peak is sim- 
ilar to the solar abundances. There is a peak around 
A = 72 and the abundance production decreases rapidly 
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Fig. 7: The number fraction relative to the solar abundances. 
As rj becomes larger, the maximum A becomes larger for rj < 
1 x 10~ 4 , while for rj > 1 x 10 -4 , the maximum A becomes 
smaller. 



above beyond the peak until A ~ 98. The maximum A 
synthesized is 98. 

To compare our results with observations such as metal 
poor star abundances, we need to take into account dy- 
namical mixing after the epoch of BBN. This depends on 
a model significantly and will be a future work. 

We should examine the idea presented in this paper 
with more realistic model, and determine whether heavy 
elements were really synthesized in BBN or not. The 
former is to take into account the diffusion effects before 
and during BBN and also lepton asymmetry. The latter 
is to calculate the nucleosynthesis in supermassive stars. 
This is because supermassive stars are generally thought 
to have synthesized first heavy elements in the universe. 
We need to know whether heavy elements observed in 
high redshift were synthesized in BBN or supermassive 
stars. Nucleosynthesis in supermassive stars and BBN 
in high baryon density region is similar. It would be 
a problem how to distinguish these two nucleosynthesis 
from observations. 



V. CONCLUSION 

We have investigated BBN in high baryon density re- 
gion. In these regions, not only light elements which 
are synthesized in standard BBN but also very heavy 
elements are produced. We found BBN is both the p- 
process like and the r-process like. The transition from 
the p-process to the r-process is due to the Coulomb bar- 
riers of proton-rich nuclei and the amounts of neutrons 
when heavy elements begin to be synthesized. The loci 
of the reaction flows change drastically above rj = 10~ 3 . 
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Above r\ = 10~ 3 , a lot of seed nuclei cause active pro- 
ton capture and the reaction flows end before very heavy 
elements are synthesized. 

Our calculations demonstrate that very heavy elements 
can be synthesized in BBN, including proton-rich nuclei. 
These nuclei will be related to the origin of the solar 
abundances, heavy elements observed in high redshifts 
and early star formations via cooling effects. 

For more realistic models in BBN, we need to include 
diffusion effects. Comparison with the nucleosynthesis 
in supermassive stars is also important. We leave these 
issues for future study. 
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